LIBRARY 

US NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 



UNITED STATES 

NAVAL POSTGRADUATE SCHOOL 




THESIS 



SIMULATION TECHNIQUES 

FOR BASIC NONLINEARITIES IN SERVOMECHANISMS 

by 

Hui S . Song 

Lieutenant, Republic of Korea Navy 



1 2ND P2238 (1-59) 



SIMULATION TECHNIQUES FOR BASIC NONLINEARITIES 



IN SERVOMECHANISMS 



Hui S. Song 



SIMULATION TECHNIQUES FOR BASIC 
NONLINEARITIES IN SERVOMECHANISMS 



by 



Hui S. Song 

// 

Lieutenant, Republic of Korea Navy 



submitted in partial fulfillment of 
the requirements for the degree of 

MASTER OF SCIENCE 
IN 

ELECTRICAL ENGINEERING 



United States Naval Postgraduate School 
Monterey, California 



19 6 0 



Sc>KlC=v k\ 



SIMULATION TECHNIQUES FOR BASIC 



NONLINEARITIES IN SERVOMECHANISMS 
By 

Hui S, Song 



This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 

ELECTRICAL ENGINEERING 
from the 

United States Naval Postgraduate School 



ABSTRACT 



Simulation techniques for nonlinearities in servomechanisms using 
the analog computer have long been subjects for investigations. In 
laboratory, these techniques play important roles for prediction of 
performance or in the study of nonlinearities in servomechanisms. 

Many successful methods in this field were developed already. 

Seven different basic nonlinearities have been chosen for investi- 
gation in this report. Net works for simulation of these nonlineari- 
ties are then generalized or initiated with complete analysis. Method 
of scaling is discussed, and -illustrative simulation in a second order 
system was included to illustrate the ease of simulation with reasonable 
accuracy. 

The writer wishes to express his appreciation for the assistance 
and encouragement given him by Professor George Julius Thaler of the 
U. S. Naval Postgraduate School in this investigation. 
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1. Introduction. 



In most existing physical systems, it is usual to have some degree 
of nonlinearity in the system, however small the degree may be. Or the 
system may deliberately be designed nonlinear such as a relay system 
utilizing the available torque to it’s full advantage* 

If this is the case, equations describing the existing physical 
system becomes nonlinear, the concept of roots of a characteristic 
equation is not defined and there Is no formal mathematical relation 
between the time domain and the frequency domain, since the principle 
of superposition does not apply, thus invalidating the basic tools for 
the analysis of linear system. For the transient response of the 
system, numbers to predict transient performance are not readily 
available. 

Thus a method of prediction of performance for nonlinear systems 
becomes very important in the field cf design of servomechanisms and 
to understand the effect of nonlinearities which may be present in 
the system; or in the design of nonlinear systems such as relay servo. 

Todate, the method of analysis of nonlinear systems is largely 
confined to phase plane method. But we are all aware of the laborious 
works and time involved with this technique as the order of differential 
equation describing the system becomes higher. And the error involved 
in the evaluation of transient performance from phase plane method in 
higher order system may become greater. 

One method of predicting system transient performance and fre- 
quency response directly at the laboratory is to simulate the problem 
to obtain solutions of equations that describe the system by use of 
analog computer. 
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Although a number of basic circuit to simulate nonlinearity have 
been introduced, a general tendency to avoid utilizing such circuits 
still remains. This report is dedicated tc encouraging the method of 
analog computer simulation in the study of the behavior of nonlinear 
mechanism which are commonly met in the usual servomechanism, by form- 
ing up the circuits, analyzing them in detail and by establishing seal 
ing. 

These circuits are formed such that it will require the minimum 
possible number of amplifiers with reasonable accuracy for engineering 
purposes. Seven basic nonlinearities are selected for the subject of 
investigation in the report, and the method of simulation is based on 
the transfer function method to minimize the number of amplifiers and 
for ease of scaling for all nonlinearities selected. Section 2 con- 
tains separated circuits to simulate basic nonlinearities with de- 
tailed mathematical analysis and general scaling. Section 3 illus- 
trates actual usage of the cirucits discussed in section 2, in a 
second order system with resultant transient 9 phase plane and fre- 
quency response curves. The method of scaling is illustrated in 
sample simulation in section 3. 
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2. Circuits for simulation of basic nonlinearities. 

In this section, basic nonlinearities encountered in the usual 
servomechanism are represented by simulation circuits for use with 
analog computer. These circuits are analyzed and simulated in the 
Donner analog computer. Statistic curves were taken and attached to 
this report. 

Scaling in general .... For scaling of a system with nonlinearity, 
it is best to neglect the presence of nonlinearity and calculate maxi- 
mum values of the desired quantities to establish a proper scaling for 
the system. In this way, the main circuit to represent the linear por- 
ticn of the system can be formed, after which the portion of the nonlinear 
circuit can be inserted. If the availability of components necessary for 
the nonlinear simulation circuit is limited, such as voltages of batteries 
or potentiometer size, it may become desirable to select scaling factors 
convenient for the formation of the nonlinear portion of the circuit. 

In that case, it would not require a great deal of time and effort to 
change scaling factors once established. 

The nonlinearity is then inserted in the linear circuit to form the 
main branch of the simulation circuit as in the case of backlash, re- 
lays and dead space etc., or added to the linear circuit to form a 
compensatory additional circuit of the system as in the case of cou- 
lomb friction. 

If the nonlinear circuit is inserted into the linear circuit to 
become a part of the main branch, the scaling of the nonlinear cir- 
cuit must be such that the scaling factors of the quantity be identical 
with that of the linear circuit at inlet and outlet stations. 
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If the nonlinear circuit is added to the main linear circuit to 



become a compensatory additional circuity scaling must be done in 
accordance with the equation describing the portion of the main branch 
to which the nonlinear circuit was added* 

These are covered in detail in section 3 of sample simulation 
of a second order nonlinear system* 
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A. 



Saturation. 

(a) Saturation phenomena in amplifiers etc.., is represented by 
input vs, output relation as shown in (Fig. l-a) 0 

(b) Circuit to represent saturation. 

i. Circuit utilizing batteries (Fig. 1-b). 



Where rdl, rd2 are diode resistances,, which further are divided: 
rf.... diode forward resistance, 
rb..,, diode backward resistance 
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The result of this circuit is the negative of the desired relation 
shown in Fig, 1-a, 

ii, Circuit utilizing power sources and potentiometers (Fig, 1-c), 



Considering the upper diode while not conducting, the cathode 
potential "P" is: 

p - f — — -t j 

- C^2.& r 6 - 40*0 

When P equals to zero, it conducts, 

A x E ( I Xd ^ o 

X, = -±- E ’ 

I “ £1 a. 

Once the upper diode conducts, 

. &z)r ^ a . 

y - — f 4r ] X| — 

Ao n i | -x ft v / i— a 

ft* 

Provided the diode forward resistance is neglegible with respect 
to other resistive parameters, 

if R,» cur ^ >> <ur 



ill « t a*r ^ | 

or R, N 1 > — 7T~ 

K' Ro 

. ^4 _ 

y a/ «. — .. £ 

~ i — a z 

And saturation voltages El, E2 will be: 

£" — ^ H- 

(c) Static characteristics. 

Static characteristic curves vere taken with the circuit of 
Fig. 1-b, where: j*. | H. 

E| - Bi. - T. 
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i. Input vs. output plot. 

This plot was obtained by use of an X-Y plotter, 
with the input sine wave frequency of .01 cps. 
and the result is shown in Fig. 1-d. The choice 
of input frequency was purely for convenience for 
the plotter. 
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Fig* 1-a Saturation characteristic 

R 0 




Fig. 1-b Simulation circuit for saturation 
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Fig. l-o Simulation circuit for saturation 
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B. Dead space. 

(a) Free play or dead space about the reference point of the load 
is represented by input vs. output relations as shown in(Fig. 2-a,). 

(b) Circuit to represent dead space. 

i. Circuit utilizing batteries (Fig. 2-b). 

let E, - E* = § = b 

MSS l) I'Ell / |E*| > X ' 

diodes do not conduct. 

Xo =o 

case 2) X\ <(~ 

lower diode conducts. 

x 0 = - &*-(- *> + e O 

provided that the diode and battery resistances are negligible 
with respect to R1 and RO. 

case 3) X| > Ei 

upper diode conducts, 

x. = - & 

provided the diode and battery resistances are negligible compared 
to R1 and RO. 

The result is the negative of the desired relation as shown in 
Fig. 2—i. 

ii. Circuit utilizing voltage sources and potentio- 



meters (Fig f 2-c), 



10 



Considering upper diode only, if the diode does not conduct, the 



cathode potential "P" is; 



- d, & -V- d, )Xj 

When P equals zero, the diode conducts. 

E -h 6~cOx, ~ £ 

*, = • E 

while the diode is conducting. 



( X' + •Cj> < p53 y . + Cx\- + C f > R, = O 



_*! g 

Xj - - 2} r h 



but 



6-^ 



"4* ’ir 

a.r Ki 



x 0 ^ - 2° 



j? s/ 

R . *' 



if 



- _ Ro 

R. 

a. ^ 

T^T« I 






Hi. -t- 2iT_ 

»-€l, «, 



«l 

At CC 



-Cx- &O 



and 



- & t 

/ A i 

pt e 

*If slope b is defined as; 



b = — 

Ri 



I 



i 



— -+ 



<3 t lr 

Ri 



Then the solution will be theoretically correct. It was assumed 
that the "a"'s were not corrected for potentiometer loading effects. 
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(c) Static characteristic curves. 



Static characteristic curves were taken with the circuit 
o£ Fig. 2-b, where; ^ = ^ w * 

- I M - 

~ i =. aid . IT. 

i. Input vs. output plot. 

This plot was obtained by use of X-Y plotter. Input 
sine wave frequency of .01 cps. were chosen for pure convenience in 
the use of the plotter. Result Is shown in Fig. 2-d. 
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Fig. 2-a Dead space characteristic 




Fig . 2-b Simulation circuit for dead space 
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C. Coulomb friction. 



(a) Coulomb friction at motor or load shaft is represented by 

the relation shown in Fig. 3-a. 

Where; XI is velocity, 

XO is coulomb force. 

The coulomb friction can be referenced to the load or motor shaft 
as desired, and the simulation is basically similar. 

(b) Circuit to represent coulomb friction. 

i. Circuit utilizing batteries (Fig. 3-b). 

Analysis of this circuit follows closely to that of the one shown 

R# 

in Fig. 1“ c except the term — - goes to infinity in this case, and the 

r\\ 

outputs are shown below, 
case 1) X\ > 0 

X 0 = " t 

case 2) X', < 0 

Xo« - E 

Where diode and battery resistances are neglected. The result is 
then, the negative of the desired relation of Fig. 3-a, 

ii. Circuit utilizing power sources and potentiometers 
(Fig. 3-c) . 

Analysis of this circuit is similar to the circuit for saturation 
except the loop gain is infinity. 

+ 

case 1) < 0 

K £ "E i if 9.?-JL at i 
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case 



2) X v >0 






(c) Static characteristics. 

Static characteristic curves were taken with circuit Fig. 3-b, 



i. Input vs. output plot. 

This plot was obtained by the use of X-y plotter, with 
input sine wave --frequency of .01 cps^, and the result i„ 
shown on Fig. 3-d. 

(d) Discussion. 

This method of simulation requires velocity to be available. 
But this is the usual condition in simulation. Due to the open feed 
back loop at the instant when XI equals to zero, the gain is infinite. 
The noise in amplifier will cause diodes to conduct one way or other, 
making output never to be zero. But this condition will not affect the 
actual simulation of coulomb friction because the physical system will 
behave in like manner and the coulomb torque applied In simulated cir- 
cuit while velocity is zero will not affect the output. 



where : 



X\ - T o iOi: 

K, ~ 1.0 M- 

e, =. Si = 4- S’ t . 
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X* 



Fig. 3-a Characteristic of Coulomb friction 




Fig. 3-b Simulation circuit for Coulomb friction 




Fig. 3-c Simulation circuit for Coulomb friction 
17 
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D. 



Simple backlash. 

(a) A simple backlash assumes two conditions. 



i. Output shaft does not drift when the gear is not in 

contact . 

ii. When the backlash is taken up and recombination occurs, 
output shaft velocity becomes instantaneously equal to 
input shaft velocity. 



These conditions automatically imply that the output shaft has no 
mass and has a great amount of friction in comparison with the inertia 
torque. The relation between input and output shaft position can be 
represented as in(Fig. 4-a.) 



XI.. Input position 



XO. .Output position 



^..amount of backlash 



(b) Circuit to represent simple backlash. 



i. Circuit utilizing dry batteries (Fig. 4-b). 



case 1) > + E 



upper diode conducts 







* t E 



y-/-* 



Provided the diode resistance is negligible. 
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y v _ -j — r - b ^ i — Yu t j 

Vlf C ' |?c sj) — j^Ts C"~bXi f 

V - ~ k*i ^ 

^ I tRCS 

Substituting eq. 4.1-4 into eq. 4.1-2 

Vi=_bX, — Yq. -f"E 

_ ~bRC S Xi j- RCS E 
I *f RCS 



V 0 « - < Yx + Y±) 

- . (*bX| j~F )Q+ ft c 5 ) 

I n- RCS 

- bX , - B ?-'-6 



case 2) ^ <- H, 

Lower diode conducts. The analysis is similar to the case 1, 
and it eventually yields; 

x 9 = -bx, -t & £/- 7 

case 3) Initially Yq_= O 

- ^ < Y, - -bx, <-» 5 



diodes do not conduct, and 



v 2 =0 



Xo - O 



After diode has conducted; 

Yq. # 0 
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and diode cuts off when; 



- 6 <- C b*t -V" Yq. ) t- £• 

1 ~ ~ (bX\ 4- Vif ) , 

\ j i US 

-VE- ^ -Clot* -t- Y4.J ' 

From eq. 4. 1-4 

* -v- Res 

or 

V 4 -+ Res Y* * - b ^ "E 

.* b*i 4-YV® E ' Re. s Y 4 . <£./— <? 

Substituting eq. 4.1-9 into eq. 4.1-8 



= R c s Viu - E 



Re S Ytj. — 0 

or - &*>-o 

/• Upper diode cuts off when Y4 reaches maximum. 

+ E = -pcsYti. + E 
Res Yq — o 

( Yv.)* O 

, Lower diode cuts off when Y4 reaches maximum in the negative 
region. 

And if R C « l 

*i — O wKe>v £L v - o 

dtb 

So, the above condition can be modified to say that when XI 
reaches maximum the diodes cut off provided RC is very small in com- 
parison with unity* When diode cuts off 




ci 

ci± 
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*3 — - constant at It’s maximum value clue to charges 

in the feed back capacitor C> and XO remains constant at the maximum 
value . 

slope b = gain constant of the first amplifier, 
backlash ^ ~ -2 H 

b 



ii. Circuit utilizing power sources and potentiometers, in- 
stead of the battery (Fig. 4-c). 
case 1) r ( > E 



Upper diode conducts 



X - 


- C lo*i *t Vh.) 


4.i~l 


4- 


y _ iL & 

— 1— . _1_ . x . 1 

CXI - ' Rt ^ 


4.2-2. 




provided /? £ OO and diode resistance 




is negligible. 




y 


y _ — e 

n i— « 

i -y _l. ajl + 

\-<V. ^ R v Rj 


4^-3 


Let 








i j. ^ AV~ <*lr 

1 + Pi -r -ft ■* ^ 




Then 






11 

>-* 






Y*-- 


" bX ‘ •+ pj-E 


■4.2-5’ 


l -+ <5~R^s 


Substituting eq. 6,2-5 in eq. 4.2-4 
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= jrC-^-Y* + jh E) 

-I?,CS fc>x, -t- R,CS 

Xo = - 'C ^ ■+■ Vq. ) 

__ P,cs^i i- f?>cS — E -b* . + Pa 
I ~+“ cP&iC S 

_ 0 ~ e) 

i -+- cf fix c s 

If <T*| (see discussion) ^ Xj> ^ “bXt "V- ( ^ 

case 2) y ( < - 1 

Lower diode conducts, and the analysis is similar to case 
and it yields: 

X. X b *i — — & 

case 3) Initially y — Q 

- E < Y, = -bx, < e 

diodes do not conduct, 

Xo - O 

Once a diode has conducted, 

Y* * o 

When a. 

-X-E: T=T“ E 

i-a 1 * 

diode cuts off. 

From eq. 4.2-5 
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4*-7 



Y = -bX, •+ Er _ 

I •+ cfftcs 

4- ^ — — bXt 4 p^ & 



-t-bX, = jppE _ 

Comparing eq. 4.2-8 and eq. 4.2-6, 
if S (?'C£ = o 

or Y4 reaches maximum value, upper diode cuts off. 

tf <r^, c « i 



<L V, = O w/h€.^ 
dt 

If Y4 reaches negative maximurp, 



lower diode cuts off. 



jj. -o '* h&K ji x ' 

slope b = gain constant of the first amplifier, 
backlash 4 = E 



provided 



s 






(c) Static characteristics. 

Static characteristic curves were taken with circuit shown 
in Fig. 4-c, where; 

X, =. txrt 

b - \ 

A® 8 ** 

C ?! m4 

- 0 *- H * 
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£ - 4!r T- 



Y - \00 K 



A = • I S-l . 



“'(this value could not be made 

lower due to non availability of 
higher voltage battery, ) 



i. Input vs. output plot. 

This curve was obtained through the use of X-Y plotter. 
Input sine wave frequency of ,01 cps, was chosen for the convenience 
of plotting, and the result is shown on Fig. 4-d. 

(d) Discussion. 

In any simulation of a' servo system, the output velocity must 
be obtainable. In this backlash case, one has to have a differentiator 
just after the output Xo in order to obtain the output velocity. To 
avoid the erroneous result of the differentiator, one usually has to 
have some small resistance connected in series with the input capa- 
citor if the input is of the step nature. Once a resistance is connect 
ed in the input side of the amplifier, the differentiated quantity will 
have error. If such is the case and the RCS in the denominator could 
be tolerated, the out put velocity is directly available in this 
simulation circuit. Y2 in Fig. 4-b and 4-c is the one, ie : 

From eq. 4. 1-5 






- b«cS X, 



7 Rcse = c? 



i -V RCS 



if 



RC «| 



Yi Z -bRCX { 



(while diode is conducting) 
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Y — O 

* Z ~ (while diodes are cut of ) 

So, the voltage Y2 represents approximately the output velocit 



The quanity 



<Ts t - & 






air 



S' ~ i if — « i 

l~ 6*" 



ajr 



« i 



« i 



(It has been assumed that the potentiometer loading effect is not 
corrected. ) 

The above condition can be realized with the use of higher vol- 
tage sources and lover resistance potentiometers with respect to the 
value of f?l. 



Also, if a corrected potentiometer setting is chosen, and if 



iL « \ 

l-O. 

S' ~ i 
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Y ► 



Fig. 4-a Characteristic of simple backlash 
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E. Ideal relay. 

(a) An ideal relay has no dead space and acts on an infinite- 
simal amount of the error signal. It can be represented by the re- 
lation shown in Fig. 5-9. 

XI. . . .Error signal. 

XO.... Relay output voltage 

This curve has exactly the same relation between input and out- 
put as coulomb friction. The only difference is in the quantity that 
input and output represent. For this reason the circuit used for simula- 
tion of the coulomb friction is also used as a simulator of an ideal re- 
lay. 

(b) For circuit; analysis of it, and static curves, refer to para- 
graph C of coulomb friction. 

(c) Discussion. 

Due to the infinite gain in the open loop, the output voltage of 
this circuit will be El or E2 even if the input is zero, which seems 
to be contradict ional to the characteristics of the ideal relay. 

In order to remedy this situation, reducing the loop gain has been 
tried. But then also, this will not be a true simulation of an ideal 
relay because of the finite slope in the output voltage to reach the 
desired level of El or E2 as illustrated in Fig. 5-6. 

If El and E2 are small, and if the gain is sufficiently large 
so that the finite slope is tolerable, there are three ways to obtain 
such a high gain. 

i. Decrease Rl, (Refer to Fig. 5-c.) 
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This method seemed to be the best way to get high gain without in- 
troducing more noises among three methods. But R1 can not be decreased 
too low if voltage sources and potentiometers are used in place of 
battery as shown in Fig, 5-d. 



i f «^»i say about 500; 

X„ = 0Tx, 1 

ay R ' 

and may become not negligible. 

ii. Increase R0. 

Since the maximum value of resistance available is 
limited, writer had to connect many resistors in series in the feed 
back path. This situation was rather an awkward one. It also tended 
to increase the noise effects. 

iii. Fractional output voltage feed back method. (Refer 
to Fig. 5-e) . 






-ih J_X| 
K, cv 



With the potentiometer coefficient "a" very small, the gain can be 
boosted to a high value. But, for some reason, this method had more 
noise effect than either of the two methods discussed previously. 

Taking these factors in consideration, plus the fact that the 
physical system will also have noise such that the output will never 
become zero, writer has chosen the open loop circuit for simulation 
of an ideal relay. 
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X\ 



Fig. 5-a Characteristic of ideal relay 




Xi 



Fig, r j-h Simulation circuit for idoal relay 




32 







. 
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Re 




Fig. 5-d Simula i, ion circuit for ideal relay 



Rc 




Fig. 5-e Simulation circuit for ideal relay 
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F. 



Relay with dead zone. 

(a) Relay with dead zone bat no hysteresis effect, has the input 



vs. output relation as shown in Fig. 6-a. 

XI..., Relay input voltage. 
XOo... Relay output voltage. 

Do . . .Dead zone. 

(b) Circuit to represent the relation of Fig. 6-a. 
i, Circuit using dry batteries. (Fig. 6-b), 



case 1) — £ < X\ < -t- E 

Xo - ° 

case 2) -£■ = ^ 

Xo= (X -§?(*•-&) ^ ^ * 

if XO becomes El; 

*. = - Cx, -H)-£, 

If JMt ~ 0 

R o + r* 

X. ^ 

T- £. 

case 3) X| <£ ~ £= ~ ~~ ~j[, 

Y, „ + E 

*■> - I; (X, -£) until -(*,-£) =* E * 
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. 














if yp ^ 



r\t 






fc hf 
Rot ff 



R 



(Xi~£) + £j. ~ 



provided -"dr O 



t± 

/?» 



Result of this circuit will be the negative of the desired 
characteristics if ___£ >■> / 

R, 



ii. Circuit using pcfwer sources and potentiometers. 



6-c) . 



case 1) 



a. 



l-O s 



E b < X, < ~E i — 



a 



D 



I — 0 : 



K- 



i 






~ C *■ - - 
0 



& s e 0 



Let 



Then 



S*t 



Rc 
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-t- 






1— <13 




R» 






-f- 


3? ^ 




<-«3 




K. 




Si, B 


5) 




a 
1 — 


l-qj 


;> (§3 


1 

* l 




2 * E s] 

!-<i* ^ 


(until Vc? 





Qj c 



C\q. 

I — Oq. 



Eq. 



If jXtj reaches — Ei± 

1-^9- 



(Fig. 
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Then 



</ 



x. = t ( x _ +■ 2i£~. e 

A,> - R, 1 .-a; 5 ' - l-tf* 

/v _ 

provided r 

case 2) _ £- ^ ^ v / r- ^3 

'^ 7^3 < *• < 

ru o Xo - o 

The result of this circuit will approximate the relations shown in 

Fig, 6-a in the negative sense. ( ^ does not affect the characteristics 

o 

of this circuit if the condition j is met.) 

R$ 

(c) Static characteristics. 

Static characteristic curves were taken with the circuit of 
Fig. 6-b, where: 

X\ — Ho 4>w wfc 
§ s S’ v* 

= £z = 7 ~ 

8 3in -**«< 

R\ - • I M 

^ of 0 



feed back pot. setting a = .5 

i. Input vs. output curve. 

This curve was obtained by using X-Y plotter with input 
sine wave frequency of .01 cps. as shown on Fig. 6-d. 

From this plot, it can be seen that the diode cut off 
and conduction voltage has a little difference. If a person wishes to 
correct this effect, it can be done by using feed back as is shown in 
para. G of relay with hysseresis. 

(d) Discussion. 

Aside from errors arising from noncritical components of the 
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circuit, the most important factor in this circuit for successful 
operation with reasonable accuracy is to have the loop gain as great as 
can be obtained. But the boosting of this gain is limited, because of 
the dead zone associated with relay. If the gain is too high, simula- 
tion shows unstability in the dead zone due to noises amplified. If 
the gain is too low, it would not simulate the desired characteristics 
of this relay because of the finite time required for output to reach 
the relay output voltage. So, there is a compromise between these two 
factors . 

The degree of compromise will depend on the computer and the 
problem itself. The gain was chosen as 400 in this case from the obser- 
vance of the performance of this circuit on the brush recorder. The 
output voltage jumped almost vertically to relay output voltage from 
zero without showing noises in the dead zone, with this gain on the 
brush recorder. The input voltage was chosen a low time rate of change 
voltage wave for this pbservance, it was .01 cps. sine wave of magnitude 
25 volts. 
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Fig. 6-a Characteristic of relay with dead zone 




Fig. 6-b Simulation circuit for relay with dead zone 




r 
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G. Relay with dead zone and hysteresis. 

(a) This relay will have unequal "pull in" and "drop out" vol- 
tages, and it has the input vs. output relation shown on Fig. 7-a. 

XI.... Relay input voltage. 

XO.... Relay output voltage. 

■^...."pull in" voltage. 

^...."drop out" voltage. 

(b) Circuit to represent relay with dead zone and hysteresis, 
i. Circuit utilizing batteries (Fig. 7-b). 



case 1) + E ' > X, ~ £ 

Starting from zero, 

Y,- 

X . = o 



case 2) 

(1) When XI increases from zero in the positive direction. 

Xi > + E 

Y, 

r z = + e r 



where voltage drop in the diode is neglected. 



Then 



x. . $<*>-!> 



Xo O 



until Xo reaches El. 



r, =- cx, = &/*» “ vbXo ^ ^ 









, -4b * 



Ri 



7 -/ 
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For the desired condition that the Xo jumps Instantaneously to 

o 

El or E2 from zero, gain — ?- .»/ as was discussed in previous para- 

Ri 

graph. If the amplifier of the computer is a perfect one with no inter 
capacitance and a perfect 180 degree phase shift, the output Xo will be- 
have as indicated by eq. 7.1-1 during the fraction of a second when 
Xo is in the jumping stage to El or E2. 

But, due to the input and output diode limiter, current can not 
reverse in that same diode, and in the limit, eq. 7.1-1 will behave 
like eq. 7.1-2. 



Xr 



- £«(*,- 

Ri 



P_ 

2 



7 J-i 



(until Xo reaches El or E2) 

But, an amplifier can not usually be perfect, then it is most 
probable that oscillation may occur during this stage, making it un- 
stable. However unstable it may be, oscillation can not occur in 
this circuit because of two sets of limiting diodes, and this condi- 
tion will tend to make a steeper slope when Xo jumps to relay output 

voltage from zero. When Xo reaches El. 

J±£o 

x 0 = •+ Bi 

v 

where diode forward resistance. 



since 



« I 




+• £r at X = 2 -Awhile XI is increasing. 



(2) After XI has reached the maximum and decreases in the posi- 
tive region. 

x.= rx, + f ) tE. « e» 

Rs 

M = Xi -t y e i 
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input diode will cut off, 




or, P n 

X, tbE.^i ^ e Xi<f-bE, 

a 

If this cut off voltage is to be 




£__d 

b~ — — 

X 0 — E tr, f or X, > when XI increases. 

Xo — + E | until X, s when XI decreases. 

case 3) X x <f- *E 

Analysis is the same as case 2 except that the other diode con- 
duct and signs reversed. 



Xo ~ ~ 2:<> for when increases. 

X* ~ -^unt il when decreases, 

ii. Circuit utilizing power sources and potentiometers, 
(Fig. 7-c). 

Starting from XI = 0 and increases in the positive region. 

X ( = o 



for - E, <* t < 



Qy 



when 



1-^3 

X z =0 
Xt> = O 
X, > 



1-4: 



■I* 



03 

I- 



n- 






3i_ + 

'-«i ic. 



fx,- liLt 

"1 i-c»j 
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Let 






a 



-t Zz* 



R e 

K. 



& v 



But, if ^ ^ oOj Xo will reach instantaneously to 



C\ ^ 






and Xo will be limited at that value. The possibility of oscillation 
due to positive feed back is eliminated by diode limiters as discussed 
in the previous pages. 

When Xo reaches instantaneously to: 



x. = 






CKa r: 
tr^ 

|-<W 



Rj> QjtK 

go TJH fy f ^ -f 
Ri 



X 

G 



. C|3 _ ^ 

F . _ — — E? ] ■+■ 

,35; * >' a s v 






Eu 



if 



TT» 



« 



« i 

K, 






Ea. 



After XI reaches the maximum and decreases in the positive region, 

<*3 

the diode in the input limiter will not cut off atX,= i-q ^3> where it 
started to conduct while XI was increasing, due to the positive feed back; 
ie., a variable DC bias voltage. And it cuts off when; 

0.1/. d i 

X|+fc PJ» E ‘ = ~s t S 



or 



Q 3 r . a <+ _ 

x,= —J*- h TT7T ^ 



If 



4* 

2 - 

, av. P _ 

This circuit will give similar results as the circuit appears 



on Fig. 7-b, provided following conditions are met. 
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(c) Static characteristics. 

Static characteristic curve was taken with circuit appear- 
ing on Fig. 7-b with input sine wave of frequency .01 cps., where; 



b - 3 different values. (setting of b was accomplished by 

setting potentiometer coefficient of Xo, while the gain 
of first amplifier was maintained at unity). 

Magnitude of input sine wave was chosen two different values to 
see the behavior of the circuit while Xo jumps from zero to El or E2, 
on the brush recorder. 



This curve was plotted by X-Y plotter and is shown on 
Fig. 7-e, where; 



g ^ I ^ TT. 

E, = Et - aa.fr r 

(?,= '! M* 

R, =. a o H- 



i. Input vs. output plot. 



^ _ qS~ A** '-0 £ 
^ - .o\ 



Fig. 7-3(1) 





Fig. 7-e (II) 





Fig. 7-3(111) 



L-<* 




2. i 



ie. a relay with dead space only. 
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(d) Discussion 

i. Since the method is based on using a variable DC bias 

voltage, the success of the circuit will largely depend on the gain of 

Ro 

the second amplifier -rr- * 

K \ 

The greater the gain the better the simulation will be. But this 
gain is limited. Because, for higher gain, the noise amplification be- 
comes great and it appears in the output signal in the dead zone, 

o 

In this case gain — % 200 was chosen and the result was satis- 

P c( 

factory. However, there was a constant error in -L _ , as indicated 

A 

on Fig. 7-e. 

To investigate this error, one more plot with b= 0 (without feed 
back) was taken as appears on Fig. 7-e(III). As can be seen from that 
figure, it also has about the same amount of error, and the writer be- 
lieves that the error is most likely caused from the nonideality of 
diodes. 

ii. The error discussed above can be corrected. Fortunate- 
ly, the error is always outward from the origin, ie., the diode conduct- 
ing point is greater than the cut off point. Once this error is known 
to be present, it can be corrected by modifying the potentiometer co- 
efficient setting "b". 

iii. A very low time rate of change voltage sine wave input 
was tested in this circuit to see the slope of Xo when it jumps to out- 
put voltage level, on the brush recorder. The result was a vertical slope 
even if the input was low time rate changing voltage. It is due to the 
unstability of this positive feed back network while the transition of Xo 
from zero to output voltage level. 
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is increased, 



Ko 

This steeper slope does not imply that the gain ^ 
and it will not give better results as a higher gain would. 

Supposing M b M is zero with the feed back loop left In place, the 
diode must cut off exactly at the point where that diode has conducted. 

But, due to the finite gain — , there will be a finite time required 

R\ 

for Xo to reach El or E2. In other words, there is a finite change re- 
quired in XI in order to let Xo reach El or E2. The steeper slope of Xo 
does not mean that the gain is actually increased. It only means that 
it shows a steeper slope inside the zone ^ as is described on Fig. 7-d, 
and it does not make the zone narrower, which is the factor that 
governs the accuracy of the circuit performance. (Refer to Fig. 7-d). 
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Fig, 7-a Character iutic of relay with rload zone and hysteresis 




Fig, 7-b Simulation circuit for relay with dead zone and hysteresis 
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Fig, 7- c Simulation circuit for rolny with dead zone and hysteresis 




Fig. 7-d Period of unstability of the simulation circuit 
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3. Sample simulation of nonlinearities in a simple second order 

system, 

A simple second order linear system of damping factor = is 
chosen to be simulated with each nonlinearity discussed in section 2. 

The method of scaling associated with each nonlinearity is illus- 
trated at each nonlinearity. 

Simulation of the linear system is also included in this report for 
comparison purposes and as a basic linear circuit to be simulated with 
nonlinearity. 

Method of simulation of tfie linear system used is the transfer 
function method although the direct component method also is adaptable. 

All scaling factors and circuit parameters in the transfer func- 
tion set up of the linear system are kept consistant throughout all 
simulation of nonlinearities. 

In order for the above stated object, ie., to keep the transfer 
function of the linear system inside the closed loop to be consistant 
throughout all nonlinearities simulation, rearrangement in the basic 
system gain constant was made in the relay servo simulation. These are 
discussed in detail at paragraph F of section 3, 

Also, a. step input only was used in this report, and for each 
nonlinearity, three different plots are obtained. They are: 

i. Transient plot of output position and output velocity. 

In this plot, two different magnitude of step inputs are used with 
two different amounts of nonlinearity, and the recording was made 
through a brush recorder. 

ii. Phase plane plot of output position vs. output velocity. 
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In this plot, two different magnitude steps were used with two 
different sizes of nonlinearity, same as ones used for plot "a", 

X-Y plotter was used for plotting. 

iii. Frequency response curve. 

Two different amplitudes of sine waves were obtained from a func- 
tion generator and the signal was fed to the simulation board. For 
each signal, two different size nonlinearities were simulated. In 
this plot, input and output transient curves were obtained through 
the brush recorder from which the ratio of output to input magnitude was 
taken and the curve was plotted on semi log scale graph. The main 
object of the frequency response curve in this report is just to show 
that the f^pu) changes in a system if the magnitude of input sine wave 
is varied due to the presence of nonlinearity in the system. Although 
a relation between the magnitude of input signal and the frequency 
response curve could be established in most cases, it is beyond the 
scope of this report and is omitted. 
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A. Linear system. 



(a) Block diagram (Fig. 1-1) 

(b) Mathematical analysis. 

£ -r & - K £ 

Oo = ^ 

Dr 5t i" s + + o;^ 

= i.r 

* - t 

CUrH+p = 6 

* 

~fc p — a2. 

Mpi = 1-33 

For step input ^ - y u(fcj 

± ^m^)=r . 5 G c 

0o - . 9 £ 7 > 

Let r £ 

Then 

t P ^qc . e^ x = i.i c m 

~^( O *&*) — - *? 2 '1 ~ • S j 

(c) Scaling and circuit (Fig. 1-2). 

This circuit will be used for all other nonlinearities simu- 
lation without changing circuit parameters. 




Oo 
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d-t $ 7 - 

And let: 

<f e. = <* e? .- , . 0 ,r 

</© = -lii = , eii 
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(-4SS) R+Cfl- = . 


5085 




sI'OO M 


C^< f ~ I- o ( huf 




= C-353j>f?f o3 Cf5 






; |.0O M. 



(d) Discussion. 

For a linear system, the frequency response must be identical 
for all sizes of input sine waves. The frequency response curves 
shown on Fig. 1-5 verifies it with little departure. This deficiency 
was believed to have arisen from reading off the values on the brush 
recorded curves. 
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Frequency response data 
(linear system) 
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Fig. 1-1 Block diagram of the basic linear system to which 
various nonlinearitios will be added. 




Fig. 1-2 Simulation circuit for the basic linear system. 
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B. Saturation. 

(a) Block Diagram. 

Saturation in the amplifier can always be referenced to the error 
signal voltage. If this is done, the circuit to represent the satura- 
tion can be placed in the error channel as shown in Fig. 2-1. 

(b) Scaling. 

Let the amplifier saturate at: 

t. £= 

u. 

The output of the saturated quantity is still position error. So, 
the saturated quantity must be scaled the same as positional error. 

e, «e,= £, 



Let : 



c / ^ -d O - ‘O 









— 3 o V~. 



' 2 -air 

(c) Simulation circuit (Fig. 2-2) 

Where : 



R. = = Rf = \'°° M * 



- 4* = 


. 3 3 3 




£ * ' 




_q, 




« * 

n 




-£*-t 


q.r, 


0 333 




a. 
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(d) Resultant curves. 

Fig. 2-3 and 2-4 shows transient responses of the system to 45 
volt and 60 volt step inputs respectively. The lower the saturation 
voltage, the lower the overshoot and the lesser the output velocity. 

These are well shown in both transient and phase plane curves. 

If the input signal increases, saturation plays its part for a 
longer period and the above mentioned effects become pronounced, as 
is shown in Fig. 2-4. 

In the frequency response curve, saturation decreases by 

comparison with the linear system and the saturation effect is greater 
for a larger input signal in identical system, thus decreasing the Mpu). 
This situation conforms to the established theory. 

(e) Discussion. 

i. Error which may have been introduced due to the insertion of 
a saturation circuit will be due to the term 9 *°3 3 3 

Roughly this will give 3% of error in error signal. This error 
can be made very small either by using smaller resistance potentiometer 
or by chosing larger voltage sources. 
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Frequency response data 
(saturation) 
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Frequency response data 
(saturation) 
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Fig. 2-1 Block diagram of the system with saturation. 
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C. Dead space or Free play. 

(a) Block diagram. 

A dead space is assumed to be present about the zero position in the 
spring loaded load shaft, and the output position is measured at the 
load shaft (Fig. 3-1). 

(b) Scaling. 

The input and output of the dead space circuit are output position of 
the system, so that the scaling of dead space circuit must be done in 
accordance with the output position. 



D - D*D 
0(0 £ (X & 

— D 

0 * ote 



Let ; 



I 

* * 
I I 



o =• -os 

d = * i 



o - 3-^V 
d - r 



(c) Simulation circuit (Fig. 3-2) 
Where; 

)r = >o o K 

R,- Ro» i- ooM 

E = 4- ^ v 

• a, = • o3S‘ r ] 

^ ?i a,* - a**** 
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C = |-o=| ^ 

f? 2 = .02. M Rj = 4-. 0 2 M 
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Since : 
But 
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___ $+ d(tc) 

0 O =- ^ 
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CX £»? 
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(d) Resultant curves. -- 

Since the dead space was assumed to be present about the zero 
position of spring loaded load shaft, it plays little effect on the 
system response for step inputs of the magnitude selected. 

The effect is almost indistinguishable on transient response curves. 

In the phase plane plot, it is better shown. If the size of dead 
zone becomes higher, output velocity at the start becomes higher. After 
some time has passed, the phase plane becomes almost same for different 
size of dead zone. 

However, it has a clear effect on the frequency response curves show- 
ing definitely different responses to different sizes of the input 
signal, thus indicating a definite presence of nonlinearity in the system. 

(e) Discussion. 

The slope of input vs. output of nonlinearity simulation circuit is: 



Ro I 

I rfe % 

as discussed on paragraph B of section 2. 
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However, the potentiometer loading was corrected in this case, then; 



slope = 
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The Vo error is about three to seven But it can be reduced by adjust- 

ing 4°so that £^2('|-4A = slope, 
ft Ri y 
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Fig. 3-1 Block diagram of the system with dead space. 




Fig. 3-2 



Simulation circuit for the system shown in Fig. 3-1 
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D c Coulomb friction,, 

(a) Block diagram (Fig. 4-1). 

The system open loop transfer function F(s) was chosen to be: 

4.1 



p(5) = 



I 



I + £ 

From this equation, the system equation without coulomb friction 



becomes : 



§ D +|e. 






If coulomb friction C is present at the motor shaft, equation 4.2 

4 3 



must be rewritten to read: 

*0 

4 - 

ZJ 



. £ • 
e o 4- ~ Qo -+* y - K' £L 



■£ d / 

Since -j « | for the chosen system* and if — A Q, ^ equation 4.3 can 

be rearranged to: 

( st 1} - K € - £ 



where initial condition was assumed to be zero. 

S j-tfvx o 

or 2* ' 

£> Qo ~ Q> 0 



From equation 4.4, it is obvious that the transfer function method of 
simulation can be used for systems with coulomb friction if proper scal- 
ing is done. Then the block diagram can be established as shown in Fig. 
4-1. If the coulomb friction is present at the load shaft instead of 

motor shaft, coulomb friction C can be referenced to the motor shaft 
, A 2 

by letting C - —3.C » where kg is gear ratio. 

J 

(b) Scaling. 
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Let the referenced coulomb friction torque C" be; 



i c '=.06 

;i d - -fj. 

Where it is assumed that the unit is in accord. 

Z'dc! = c' 

where 

One can choose any convenient value of C'f or more convenient values of 
oCC # Usually C' is chosen for a convenient value of limiting voltage 
of the coulomb simulation circuit. Batteries may well best be used in 
the limiter, because of the scaling convenience. 

Going back to equation 4.4 in the problem form, 

K£ ~ cVf 1 S o 

£ ■»- i 

This equation can be written in terms of the computer values as equation 





4.5. 






£ ’ c(t 

where, in this example: ^ c{ & 

From equation 4.5, the scaling procedure is clear and self explanatory, 
(c) Simulation circuit (Fig. 4-2) 
where : 



Let : 



Y - loo K 
Rx'kr 2 M - 
[?* * i. ooM 

£ - \ 00 V 
o( C ' . ol 2. 
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c’ 


c' - .06 
C'- •'! 




5- - 


II 


C 7 = 


/ 0 - - 








t 


a, [ 


~ Ft? " 






H 

a& | 
ft i ' m&x • 


0, - _ 
Erie “ 

= .oo^4 


• of 





4, Z r 

i—rti 



l-(h 



Rearranging equation 4.5, 

£ - Jw- L * - c_Vc ^ J 

Let 

cu c *~. = ./?/<? 



*7 



6d4~ 

Also, 

b_ 

° 

(d) Resultant curves. 



~t~rr ~ b- .I8 r ii f? a * »•« ° M * 



The effect of coulomb friction in the transient response of the 
system is as if the damping of system is increased. In Fig. 4-3 and 4-4, 
the overshoot becomes less when the coulomb friction is increased, while 
at the same time ouput velocity decreases. 

Due to coulomb friction in the system, the output position will have 
a finite steady state error when step input is imposed. This steady state 
error will be; 0 < £ <<s < JL! 

5:5 N K 

S = S' ^ 

~c'- 10 £ 5S $ 4<<4r v~. 

This effect is well shown in the phase plane plot of Fig. 4-5. 

The frequency response curve shows a definite presence of nonlinear- 
ity in the system. 

(e) Discussion. 

In the case^ , and the output of the limited quantity 

has a error of *£ «• 00 4- Q 0 
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This error could be reduced by the use of dry batteries in place of vol 
tage sources and potentiometers* Also the same objective can be accomp 
lished if lower resistance potentiometer or higher voltage sources were 
used* 
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Frequency response data 
(coulomb friction) 



C'~ $ 



U) 


0, 






M 


.0943 


19.2 


.96 


39.5 


.988 


.1885 


20.2 


1.01 


41 


1.025 


.283 


21.2 


1.06 


44 


1.1 


.377 


23.7 


1.185 


48 


1.2 


.472 


26 


1.3 


54 


1.35 


.566 


28 


1.4 


59 


1.475 


.613 


28.5 


1.425 


60 


1.5 


.66 


28 


1.4 


59 


1.475 


.754 


25.8 


1.29 


55 


1.375 


.848 


21.2 


1.06 


45 


1.125 


.943 


16.5 


.825 


35.5 


.888 


1.885 


3 


.15 


6.5 


.1625 


3.77 


.7 


.035 


1.4 


.035 


5.66 


.26 


.013 


.6 


.015 


7.54 


.14 


.007 


.3 


.0075 


9.43 


.05 


.0025 


.14 


.0035 






U)£* 




q. o u) 



86 



4 

01 

02 

03 

04 

05 

06 

065 

07 

08 

09 

1 

2 

4 

6 

8 

0 



Frequency response data 
(coulomb friction) 
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Fig. 4-1 Block diagram of the system with Coulomb friction. 




Fig. 4-2 Simulation circuit for the system shown in Fig. 4-1 • 
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E, Simple backlasho 



(a) Block diagram (Fig, 5-1) 

Backlash is assumed to be present in the load shaft of the system, 
and the output position is measured at the load shaft. 

The simple backlash implies that the load inertia is negligible 
compared to the drag force so that the load does not drift during 
backlash period and when gears recombine, load shaft velocity jumps to 
that of the motor. 

(b) Scaling, 

The scaling of this problem is really simple. Since the quantity 
to be dealt is the position, and the circuit is to be inserted to the 
main branch, the backlash A must be scaled to the computer value by 





l| A ~ • 



(c) Simulation circuit (Fig. 5-2). 



Where ; 



& - is- 
R — •! M. 



\ M-P 




= -035-7 



a K = .06% 8 



2, 



From paragraph D of Section 2 
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a- 



OfR c S)(~O l i‘ 1-4, 

I "f* <S RCS 



S £ I -h A + a.h , Gj ^ 

1 '-«• R7 + ff a 



Since potentiometer loading effect due to R1 and R3 are taken care of; 



S ~ — -a. N 0*4- 

i A i 



a, = ' * 3 b 'I 

a, - • 0<; ^ 



- i. ons- 

&.= fl, + ^-E 

l-o i 

f? x = .oiS- M Cj, =■ |. 0 3 4- 

(9 ^ 

b« **** -. 3 <? r 

(d) Resultant curves. 

It is a well established fact that the backlash has a destabiliz- 
ing effect on the system performance* From Fig* 4-3 and 4-4 of tran- 
sient response curves, it is seen that the overshoot and output vel- 
ocity becomes greater with more oscillations when the amount of back- 
lash is increased* 

Fig* 5-5 and 5-6 of phase plane plots, show more clearly of this 
situation* For the system damping factor of *333, of second order 
system, there should be no limit cycle, whereas the simulation shows a 
limit cycle for reasons that are not quite clear. 

However, with the understanding of the relation between maximum 
damping factor that will give limit cycle in the system, this circuit 
can be. utilized for general engineering purposes* 

Frequency response plot on Fig* 5-7 and 5-8, shows that the in- 
crease of backlash corresponds to a greater Mpu) * Mfu) *- s generally 
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higher than that of the linear system which is shown in Fig 0 1-5. 

For larger input signal, backlash plays less effect on the system, 
thus becoming close to linear system performance* which in turn will in- 
dicate a decrease in These effects are well shown on the frequency 

response curves. 

(e) Discussion. 

io If E in the Fig. 5-2 is increased to higher value, and if the 
potentiometer loading has been corrected, & can become closer to the 
unity. 

ii 0 As resultant response curves indicate, the stimulation shows 
limit cycle contradicting to the behavior of the physical system. 
Although the reasons are not quite clear, the writer suspects it to be 
caused from stored energy in the capacitor C of Fig. 5-2. 
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Frequency response data 
(simple backlash) 
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Frequency response data 
(simple backlash) 
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Fig. 5-1 Block diagram of the system with simple backlash. 




Fig. 5-2 Simulation circuit for the system shown in Fig. 5-1. 
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F. Ideal relay. 

(a) Block diagram. 

The ideal relay is placed in the error channel as shown in Fig. 

6 - 1 . 

In all previous cases, motor gain was assumed to be included in 
the constant K — .2,2 Is „ 

But, for simulations of a relay servo, it is assumed that K = 2 2 S’ 
does not include the motor gain constant, and it describes the open loop 
gain of the system not including motor. Also the motor is assumed to be 
an ideal torque motor. 



The equations describing the 


system become: 




JS 0 W 6> 0 = +■ T 


e > o 


6* 1 


j e 0 * -f Oo = - t 


£ < C? 


<$• 2 


where T denotes the constant 


torque output of the motor. 




-f s J- oTIr 






Q + 6 = 2. a. ir (t) 


£ > o 




e +6 = Z'Z-S-f-T) 


1 < o 





Then the transfer function set up used for simulation of all cases 
in previous paragraphs can be used without changing scaling factors and 
circuit parameters. 

(b) Scaling. 

Assuming that the torque T is in proper units the system equation 
can be arranged as follows: 
i. £ > O 

s & 0 cii-s) = - 2 * iS ' t 
G 0 = _^1±U 

t £ 
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£ <o 



ii , 



©. -- — ~ (- T) 

I i- S 



And the computer equations are; 

= — iAf-C4X)f ^-2 

S+^. <* + ' 

From these equations, it is seen that the scaling is very simple, 



similar to the coulomb friction case. 



If 



Ot T - o( 0 



= . oi$- 






linear system simulation circuit can be used without rescaling. 



Let 



T 



• * 
1 1 



t = • 4- Jr 
T = .qs- 



T - 30 V. 

T" =■ £To IT 



(c) Simulation circuit (Fig. 6-2) 
where ; 

R, = M - 

f = i o o K 

jr - XOo T. 

a ' = * 

" r'x b - t « So 

I ^ V ' 



|?v 

(d) Resultant curves. 



A/ 



= .13^ 

Q, * * 2 



Theoretically, the system must be zeroed in to the origin in the 
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phase plane, The simulated results shows limiting cycle contradictory to 
the theory. This is due to noise effect in the amplifier since the loop 
gain is infinity. 

In terms of transient response, the stability is poor. The stability 
becomes poorer if the motor torque is increased, exhibiting higher over- 
shoot and output velocity. These effects are shown in both transient 
response curves and phase plane plots. 

The frequency response curves are ploted as shown in Fig, 6-8 and 
6-9. The portion up to the resonance peak in the curve, is not a true 
response, because the output was an oscillatory wave. The nonlinearity 
is well verified to be present in the system by noting different res- 
ponses of the system to different sizes of input, 

(e) Discussion, 

Although the second order system with damping chosen must not 
exhibit limit cycle theoretically, the physical system will have noise 
involved, which in turn will drive the motor even when the error signal 
is zero. Considering this fact, the limit cycle observed in the simulat- 
ed circuit is not a detrimental factor. If a resistor is shunted across 
the feed back path, it eventually zeroes in to the origin, but not with- 
out a finite time sacrifice for the output to get up to relay output 
voltage. The writer preferred the infinite gain open loop circuit. 
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Frequency response data 
(ideal relay) 
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Frequency response data 
(ideal relay) 
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Fig. 6-1 Block diagram of the ays tom with ideal relay. 
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Go Relay with dead zone* 



(a) Block diagram (Fig, 7-1) 

The gain of amplifier Ka has only one effect on the system. It 
determines the sensitivity of relay dead space in terms of actual 
system position error, Thus*, if the relay dead space is referred to 

a dead space in terms of error*, Ka can be omitted. This implies that 

if the dead space simulation cirucit is placed before the amplifier 
with a proper scaling*, the amplifier to simulate Ka can be omitted. 
Let p £ width of dead space of relay in relay voltage. 

D width of dead space of relay in terms of position. 



T -■^5' 




The system equations become 



rf 6 - t 
J© -T 

J e -t -f 6 =■ o 




(b) Scaling. 



=f = 3 O TT. 



Since D is in terms of error 



; o' =-9 o— D - V~ 

\\ D' = l.? 0- * ^ O = S V 



(c) Simulation circuit (Fig. 7-2). 



Where = |, o o M 



f = \ o o K 
E, = H-o v. 
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£ z — 2. oo v 
Ro & 10M 
a « a r 



fey* o = 3 
ioY - 6 



3 



(d) Resultant curves. 

With dead space in the relay, system becomes stable. Wider the 
dead zone, the more stable system is obtained. These effects are 
shown in Fig. 7-3 and 7-4 of transientresponse curves. Because of 
noise amplified in amplifier, the limiter loop gain could not be gotten 
higher than 65 which is a little too small. Even with this gain, noises 
showed during the transition of dead space. Also in transient curves, 
drift of output position due to the noise effect is observed. 

In the phase plane plot of Fig. 7-5, two plots for different size 
of nonlinearity seems to fall on top of each other. This is due to 



gain of limiter amplifier 



K £ S’ 



i — bi 



i — 



= T 



b,~ 



btl" 



2i£ *.o I 



; J2. = 3 -r 

I 2. ^ 

;; jc_ . 6 _ 

2. 2 ~ 



I- A» 






1^4, 



a. - - o 3 <6 < 3 









b,lr 



|i- 



'R 






o ^ 



^ . o I 



^ ■ ■ 1 
0 ~ I _ 



b, E. 



b, 



3 < £ <■ 
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similar Isoclines of phase plane before and after the dead zone. As 
can be seen, two plots depart perceptively after a few cycles. 



In frequency response curves, the first part of the curve up to 
the resonance peak is not a true one because the output did not reach 
the steady state condition, 

(e) Discussion, 

Due to not high enough loop gain of 65, the simulation circuit 
imposes additional error in t he simulation. For error voltages less 
than region, the torque applied is not ^ - 3 ^ but 

£5- . This error* can be reduced if the gain can be increased. 
If batteries or low noise voltage sources were used* gain may be gotten 

L V* 

higher. The term imposes another error in the circuit. If 

fZt u ._ 






the simulation will have 1.8% error in the torque. 



But this error factor can be reduced with the use of higher voltage 
sources or lower resistance potentiometer. 
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Frequency response data 
(relay with dead zone) 
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Frequency response data 
(relay with dead zone) 
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Fig. 7-1 Block diagram of a relay servo system with dead zone. 
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Ho Relay with dead space and hysteresis. 

(a) Block diagram (Fig. 8-1). 

Let; £ £ Relay pull in voltage. 

d 1 £ relay drop out voltage. 

" 5U pi 

J_ ^ referred to error. 



d A SL referred to error. 

2- 

The amplifier gain Ka has same effect as was discussed in the relay 
with dead space s and is treated in same manner. Then system equations 
become : 



3 e n- f 6 - 0 


6<i 


|£j increase in 


+ 


jetfe-T 


<■>4 


(£| increase in 


-h 


J 0 + -f e = ° 


e<i 


i£l decrease in 




J0 -V-f 6 - ~T 


6 <1 


1 £l increase in 


— 


j e i-f © - 0 


£ >d. 

c. ^ 


l£)decrease in 


— 



Let 



i. 



ii. 



T = ' 67 
3 O 



£ = 1.1 Sr r . 
, .S T. 



£ = m r 7, 
i'-ar. 



. oqr TT. 

d _ . 0 4-S- V“o 

*X 

£. ^ .c\S v. 

i = -0>lr 



(b) Scaling 

T ® 

V 44 volt dry battery was available. 
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l ' £ - £ ^ a 4- r 

a * s! f;. a 3T. 

3. dfc> 

-E. d a 3 

* i 3 

il. _ 

i ® J2 — - s T 

» i diJ ° 

3[ „ d X — | TT 

a" IdS " 1 v ' 

2 -SL ^ H. V 

-t TT T ' 

(c) Simulation circuit (Fig. 8-2). 
where; *- 



ft £ * I M 

(?0 ^ 2.0 H 





E, - i«ot 




i. 


jp = . O 


J2 _ d =• 2- V. 
-2- 




JE - Si - *f S-i. 
2- 2. 






J} ^ JL - . oU-s-U. 

T 




ii. 


f=.°9l felr 


1 - f *->■ r. 




~ = . o 9 » 

limiter loop gain = 


IL° ‘X, 2.0 0 
* \ 


(d) 


Resultant curves. 




The 


case of -£ - - 0 


; ie. a relay with dead 



included in this paragraph for comparison purposes, because T is 
different from the one used in paragraph F of this section. 

In the transient response curves , comparing Fig. 8-3 and 8-4 to 
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Fig. 8-5 (relay with dead zone only), it is evident, the relay with 

hysteresis has more oscillation with higher overshoot. This implies 

that the effect of hysteresis in relay servo is a destabilization. For 

cases of hysteresis chosen, it exhibits a definite limit cycle, while if 
p A _ 

-gr - — ~0, system is stable. 

Fig. 8-6 shows three cases of nonlinearities. If there is no hystere- 
sis effect in the relay, it zeros in to the origin as shown by curve 
"a", while others reach limit cycle. 

The first part of frequency response curves up to the resonance peak 
is an approximation. The. output did not reach the steady state in that 
part of the frequency. The purpose of these curves are to show the 
presence of the nonlinearity in the system. 

(e) Discussion. 

f?0 

As was stated in Section 2, the limiter gain ^ must be great in 
order for a satisfactory operation of this circuit. With rectified D-C 
power source in place of E2, gain could not be gotten higher than about 
50. It was then necessary to replace E2 with dry batteries to rid it 
of noises. Dry battery available in the laboratory was 44 volt. For 
this reason, it was found necessary to change the torque to T~ 
volt from previous case of relay servo simulation. 

With the use of dry batteries, the gain could be increased to about 
200 which gave satisfactory results. 

In all phase plane plots, errors in n drop out** and "pull in ,v 
voltages is 127. in this simulation determined from the phase plane plots 
assuming that the plotter involves no error in itself. This error, is 
believed to be caused from inaccurate representation of and the non- 



133 



ideal characteristics of diodes used. 



Also, the loop gain of non-infinity introduces error. In general, 
if the power sources are noise free so that the loop gain can be in- 
creased to a higher value, and if accurate representation of "pull in" 
voltages can be accomplished, it could be reduced. 

Or if the amount of error involved with the circuit is knowjlj then 
the error can be corrected by adjusting the potentiometer setting "f". 
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Frequency response data 
(relay with dead zone, hysteresis) 
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Fig. 8-1 Block diagram of a relay servo system with hysteresis 
effects 'in the relay. 
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Fig. 8-2 Simulation circuit for the system shown In Fig. 8-1, 
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4. Conclusion, 



Simulation technique of basic nonlinearity encountered in usual 
servomechanism is generalized and analyzed in detail, with illustra- 
tive sample simulation of second order nonlinear system 0 Utilization 
of this technique is apparently easier and time saving with reasonable 
accuracy for most engineering purposes,, The scaling technique was 
illustrated. The error in general is less than 5% more added to the 
error involved in the simulation of linear system, except cases of relay 
with dead zone and hysteresis and simple backlash. But those exception- 
al cases can be made to have errors within the prementioned value. 

In all cases, error can be reduced by using more ideal components 
in simulation such as batteries, idealized diodes with more complicated 
setup or non-ideal diodes, or potentiometers well balanced. It can 
also be minimized by more fortuitous choice of scaling factors as 
was discussed in each paragraph of section 3, such as making scaling 
factor of the time greater so that an error resulting from noncritical 
component of the circuit may become negligible in terms of the physical 
system values. 

The use of this technique in the study of the behavior of non- 
linearity in servomechanisms saves a great deal of time and labor 
especially if one deals with the transient performance and frequency 
response of a system. 

Given numbers for equations describing a physical system and 
numbers concerning nonlinearities, it is easier to set up circuit in the 
laboratory and simulate them to study or predict the system performance 
than drawing phase plane plots and tabulating data to further obtain 
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the transient response. When it comes to frequency response, the only 
mathematical approach "describing function method" does not give better 
accuracy, if as good. 

For reasons stated above, the writer believes that the simulation 
method in the study of nonlinear behavior is preferable to any other 
method, and this method concerning basic nonlinearities commonly met 
are studied and analyzed with Illustrative scaling in this report. 
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